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A quasi-linear theory is presented for how randomly forced, barotropic velocity fluc- 
tuations cause an exponentially-growing, large-scale (mean) magnetic dynamo in the 
presence of a uniform shear flow, U = SxEy. It is a "kinematic" theory for the growth 
of the mean magnetic energy from a small initial seed, neglecting the saturation effects 
of the Lorentz force. The quasi-linear approximation is most broadly justifiable by its 
correspondence with computational solutions of nonlinear magnetohydrodynamics, and 
it is rigorously derived in the limit of large resistivity, 77 — >■ (X). Dynamo action occurs 
even without mean helicity in the forcing or flow, but random helicity variance is then 
essential. In a sufhciently large domain and with small wavenumber kz in the direction 
perpendicular to the mean shearing plane, a positive exponential growth rate 7 can occur 
for arbitrary values of 77, the viscosity v, and the random- forcing correlation time tf and 
phase angle ^/ in the shearing plane. The value of 7 is independent of the domain size. 
The shear dynamo is "fast", with finite 7 > in the limit of 77 — > 0. Averaged over the 
random forcing ensemble, the ensemble-mean magnetic field grows more slowly, if at all, 
compared to the r.m.s. field (magnetic energy). In the limit of small Reynolds numbers 
{rj, V — >■ 00), the dynamo behavior is related to the well-known alpha-omega ansatz when 
the forcing is steady (tf 00) and to the "incoherent" alpha-omega ansatz when the 
forcing is purely fluctuating. 



1. Introduction 

This paper presents a theory that yields exponential growth of the horizontally-averaged 
magnetic field {i.e., a large-scale dynamo) in the presence of a time- mean horizontal 
shear flow and a randomly fluctuating, 3D, barotropic force (i.e., with spatial vari- 
ations only within the mean shearing plane) in incompressible magnetohydrodynam- 
ics (MHD). This configuration provides perhaps the simplest paradigm for a dynamo 
without special assumptions about the domain geometry or forcing {e.g., without mean 
kinetic helicity). We ca ll it the elemental shear dynamo (ESD). There is a long his- 
tory of dynamo theory (lMoffatdll978l : iKrause fc RadleilllQSOt [Roberts fc Sowardlll992 : 



[Brandenburg fc S ubramaman^'2005^. but much of it is comprised of ad hoc closure ansatz 
{i.e., not derived from fundamental principles and devised for the intended behavior of 
the solutions) for how fiuctuating velocity and magnetic fields act through the mean elec- 
tromotive force curl to amplify the large-scale magnetic field. Here the horizontal-mean 
magnetic field equation is derived within the "quasi-linear" dynamical approximations 
of randomly forced linear shearing wave s and flow-induced magnetic fluctuations. 

In the standard ansatz ( MoffattI 1978f l , the mean-field equation in dynamo theory has 



the functional form of 

dtB = L B + D -.WB + ... , (1.1) 



f Email address for correspondence: jcm(S)atmos. ucla.edu 
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where the over-bar indicates some suitably defined average; B is the mean magnetic field; 
and L and D are second- and third-order tensor operators (often denoted by a and /3) 
that express the statistical effects of the velocity field v through the curl of the mean 
electromotive force, V x (d x B). The dots encompass possible higher-order derivatives 
of B (which would be relatively small if there were a spatial scale separation between the 
mean field and the fiuctuations) and resistive diffusion. If v itself is steady in time, then 
(jl.ip is an exact form for the electromotive effect, and the kinematic dynamo problem 
can be viewed as an eigenvalue problem for the exponential growth rate 7 given v; in 
this case, however, there will be no scale separation between v and B, and 7 may not 
be positive. An important weakness in such an ansatz is the lack of justification for 
particular forms of L and D in time-dependent flows. We will see that the ESD theory 
provides a clear justification, and it mostly does not fit within the ansatz (jl.ip because 
the tensors are time-integral operators except in particular limits (Sec. [5]). 

The ESD problem specifies a steady flow with uniform shear S, a small initial seed 
amplitude and vertical wavenumber for the mean magnetic field, and a particular 
horizontal wavenumber k^f and correlation time tf for the random force. It defines an 
ensemble of random-force time series that each gives rise to a statistically stationary 
velocity field, and the induced dynamo behavior is assessed over long integration times 
with further ensemble averaging. 

This paper takes a general parametric view of the ESD derivation and solutions. A par- 
allel report utilizing a minimal proof-of-concept derivation for the tr eble limit of small 



kinetic and magnetic Reynolds numbers and weak mean shear is in iHeinemann et al 



(|2011q( ): the relation between the two papers is described in Sec. 15.21 Th e experimen- 



tal bas is for developing the ESD theory is the 3D MHD simulations in lYousef et 



([200830). They show a large-scale dynamo in a uniform shear flow with a random. 



small-scale force at intermediate kinetic and magnetic Reynolds numbers. Their dynamo 
growth rate is not affected by a background rotation, even Keplerian. Additionally, new 
2"'"D simulations — a barotropic velocity with spatial variations only within the mean 
shearing plane (x, y) and a magnetic field with (x, y) variations plus a single wavenumber 
kz in the vertical direction z per pendicular to the plane — also manifests a large-scale 



dynamo (IHeinemann et aLll2011M) . Furthermore, within this 2+D model, successive levels 
of truncation of Fourier modes in the shearing-plane wavenumber demonstrate that its 
dynamo behavior persists even into the quasi-linear situation for which the mean-field 
theory is derived here. Thus, the dynamo solutions of the ESD theory are a valid expla- 
nation for computational dynamo behavior well beyond the asymptotic limit of vanishing 
magnetic Reynolds number. 

From general MHD for fluctuations in a shear flow (Sec. [2]), a quasi-linear model is 
developed for shearing waves (Sec. [3]) and for induced magnetic fluctuations and the 
horizontal-mean magnetic fleld evolution equation with dynamo solutions (Sec.|4]). Ana- 
lytic expressions for the dynamo growth rate 7 are derived in Sec. [5] for several parameter 
limits, and general parameter dependences are surveyed in Sec.[Sl Section [7] summarizes 
the results and anticipates future generalizations and tests. 



2. Governing Equations 

The equations of incompressible MHD are the Navier-Stokes equation for velocity v, 

dtv + v-\/v = --Wp + B-\/B + j^W'^v + f, (2.1) 
P 
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where / is a prescribed forcing function, density p is constant, and pressure p is deter- 
mined by the constraint, 

V-v^O, (2.2) 
and the magnetic induction equation for B (in vefocity units), 

dtB + v - VB = B XIv + riV^B , (2.3) 

with 

V-B = 0. (2.4) 

An exact, conservative solution to the above equations is given by an unmagnetized, 
uniform shear flow of the form 

v^Sxey, B = 0, (2.5) 

where the shear rate 5 is a constant in space and time and e denotes a unit vector. To 
study the dynamics of fluctuations on top of the background shear flow (j2.5p . we rewrite 
the equations of motion in terms of the vefocity fluctuations u defined through 

V — SxBy + u . (2.6) 
Assume that the volume average of u is zero. Substituting (|2.6|) into (|2.ip and (|2.3p yields 

Vu + uyu + Su^By = -Vp + B ■ VB + vV'^u (2.7) 

and 

VB + u- VJB ^ B ■Vu + SB^^Ey + ryV^B , (2.8) 

where 

V = dt + Sxdy . (2.9) 

The only explicit coordinate dependence in (j2.7p and (|2.8p arises through the differ- 
ential operator (|2.9p . which contains the cross-stream coordinate x. This means that we 
can trade the explicit z-dependence for an explicit time dependence by a transformation 
to a shearing-coordinate frame, defined by 

x' = x, y' ^y- Stx, z' = z, t' ^t. (2.10) 

Partial derivatives with respect to primed and unprimed coordinates are related by 

d^, = a, + Stdy , dy, = dy , d,, = , dv =dt + Sxdy - V , (2.11) 

which shows that the explicit spatial dependence is indeed eliminated in the shearing 
frame. Therefore in shearing coordinates there are spatially periodic solutions, in partic- 
ular a Fourier amplitude and phase factor, expressed alternatively as 

X(a;,y, z,i) = Re|x(Oexp \kx{t)x + ikyy + ik^z | 

= Ylc^x{t')Gyip ik^Qx' + ikyy' + \k^z' |, (2-12) 

where the transverse wavenumber ky and the spanwise wave number kz are constant in 
both coordinate frames, but the streamwise wavenumber k^ varies in time according to 
kx{t) — kxQ — Skyt. For an observer in the unprimed ("laboratory") coordinate system, a 
disturbance that varies along the streamwise direction stretches out as a result of being 
differentially advected by the background shear flow; for an observer in the shearing 
frame the Fourier phase has fixed wavenumbers {k^o, kyo, kz). 
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3. Dynamics 

3.1. Simplifications 

Guided by the experime ntal demonstrations of the shear dynamo ( Yousef et al. 2008Ql f6: 
iHeinemann et al. 20116 ). we make the foUowing simphfying assumptions: 

(a) The magnetic field strength is sufficiently small so that there is no back reaction 
onto the flow. In this so-called kinematic regime, we drop the Lorentz force. 

(6) The 3D forcing is restricted to two-dimensional spatial variations in the horizontal 
{x,y) plane {i.e., barotropic flow with dzU = dzp = 0). (With this assumption it makes 
no difference whether the system is rotating around the axis or has a stable density 
stratification aligned with e^. For these dynamical influences to matter, u has to have 3D 
spatial dependence.) In this case the dynamics reduce to forced 2D advection-diffusion 
equations for the vertical velocity, Uz, and the vertical vorticity, Wz = • (Vj_ x u±); 
viz., 

Vuz + u±_ ■ V±Uz = vV^Uz + fz 

VlOz + M± • V^iOz = vV\i^z + ez • (Vi X /^) . (3.1) 

We use a notation for a horizontal vector as 

a^= a^e^^ ayEy . (3.2) 

Because u has no z dependence, the non-divergence condition reduces to Vx • = 0, 
and we introduce a streamfunction $ for the horizontal velocity and its associated vertical 
vorticity: 

Mi = ezXV±$, Wz=Vi$. (3.3) 
(c) Fluctuation advection is neglected in (j3.1l) . so the vertical momentum and vorticity 
balances are linear. 

VUz = I^V^Mz + /z 

VijJz = v^l^z + ez • (Vi X /^) . (3.4) 

3.2. Conservative Shearing Waves 
For linearized conservative dynamics (/ = 0, v = Q), (|3.4p is 

Vuz = Vujz = . (3.5) 

The Fourier mode solutions are 

Uz = Re { uzo e'^ } 

c^z - Re { Lbzoe'^ ] , (3.6) 

with a phase function that can be alternatively expressed in shearing or laboratory co- 
ordinates as 

(j) = k'^x' + ikyy' = kx{t)x + ikyoy . (3.7) 

The constants k'^ = kxo, k'y — kyo, Uzo, and cJzO are set by the initial conditions, and 
a tilting a;-wavenumber is defined by kx{t) — kxo — Skyot. From (|3.3p the associated 
horizontal velocity is 

u^ = ^^^^^Re{icozoe'^} , (3.8) 

where k^ — k^ + k^Q . Notice that u± (t) grows when k^ (t) / kyo > by extracting kinetic 
energy from the mean shear (an up-shear phase tilt), and it decays when kx{t)/kyQ < 
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(down-shear). As i — cxi, m_l for any fco- This shearmg wave behavior is sometimes 
called the Orr effect. 



3.3. Single-Mode Forcing 

In a quasi-linear theory the random fluctuations can be Fourier decomposed into horizon- 
tal wavenumbers, and the resulting velocity and magnetic fields summed over wavenum- 
ber. It suffices to examine a single wavenumber forcing to demonstrate the ESD process 
(c/., (|4.26|) ). When f{x, y) is restricted to a single horizontal wavenumber in the labora- 
tory frame fe^/, we have 

/ = Re { f{t)e'^' } , (3.9) 

where the Fourier coefficient / is specified from either a random process. The spatial 
phase of the forcing is fixed in laboratory coordinates: 

4>f = kxfX + kyfy . (3.10) 

The non-divergence condition on the Fourier coefficient in (|3.9p is k±f ■ j ^ = 0; hence 
we can write 

/± = /-ie_L/ , with e_L/ = ^ (3.11) 

the unit vector perpendicular to the forcing wavevector. Here k±f = \kj_f\. The forcing 
coefficient is thus 

f = ke±f + fze,. (3.12) 
Taking the cross product oi k±f with / yields 

k±fXf = k^f{he,-f,e^f). (3.13) 

This is used to define two further relations. The forcing coefficient for vertical vorticity 
is 

0^ = -ik^if X f = ikj_ff_i. (3.14) 



The spatially-averaged forcing helicity (defined byi?=^/-Vx/y where brackets 
denote an average in the indicated superscript coordinate) associated with a single Fourier 
mode is defined by 

Hit) = i Re[/* • (ik^f X /)] = Rc[/;6,] , (3.15) 

which is a real number. The asterisk denotes a complex conjugate, and we now incorporate 
a caret symbol in H{t) to be consistent with other forcing amplitudes. 

The Fourier mode coefficients fz {t) and dz (t) are complex random time series that are 
mutually independent between their real and imaginary parts and between each other, 
and they have zero means. We consider an ensemble of many realizations for these time 
series. (We will also analyze solutions with steady forcing (i.e., with / fixed in time with 
values taken from the same random distribution).) For a given realization, we generate 
the forcing coefficients from an Ornstein-Uhlenbeck processes with a finite correlation 



6 



J. C. McWilliams 



time, tf. Thus, 



ol{ti)d,{t2) 



= cxp 
= Oz exp 



£ 



-\tl-t2\/tf 
-\tl-t2\/tf 

= 0, 



(3.16) 



where £ is the expectation value averaged over fluctuations and and Oz are positive 
forcing variances. In particular, the helicity has zero mean, £ H{t) ~ 0. 



3.4. Stochastic, Viscous Shearing Waves 

We assume single-mode forcing. For simplicity we assume that the fluid is at rest at t 
The resulting solutions to p.3p - (|3.4l) are 



0. 



Uz{x,y,t) 
^z{x,y,t) 
u±{x,y,t) 



d/iG,(t,M)Re{/.(M) e'^(^)} 
dfiG,{t,^i)Re{dz{fi) e'^('^)} 

-BzX k±{t- fl) 



d/i G^{t, fi) 



fc2(t_^) 



^ Rc{i64^)e'^(^)} , (3.17) 



which can be verified by substitution into the dynamical equations. The wavevector is 
k±{t) = {kx{t), kyf) with kx{t) = kxf — Skyft and fc^(i) = k1{t)+kyf. The phase function 
(j) represents continuous forcing at the single, laboratory- frame wavenumber fc / , and its 
evolving shear tilting is expressed in kx{t). We can write it in either the sheared or 
laboratory coordinate frame: 



{x',y',t';^) ^ (kxf + Skyffi)x' + kyfy' 
(l){x,y,t;fi) = kxit ~ fi)x + kyfy = fc^(t-^)-x. 



(3.18) 



where kx{t — ^jl) 
decay factor, 



Skyf{t — /i). The viscous damping effect is expressed by the 



G'^(t, /x) = exp - / dpk\{p-^L) 



cxp 



dcfci(c) 



(3.19) 



which is a Green's function for p.ip . For compactness we can write this as an equivalent 
function of a single time difference, G,j(t — /i). 

In the first line of p.lSp . is expressed in shearing coordinates {x',y',t'); note that 
the phase of the shearing wave is independent of t' , but it does depend on the forcing at 
the time /i when the wave was spawned. The second line is the equivalent expression in 
laboratory coordinates (x,y,t). For compactness we write this below as (/)(/i), with the 
other space-time dependences implicit. 

If V = Q (hence G^, = 1) and the forcing is applied only at the initial instant {i.e., 
fz = S{p)uzo and = S{^)(2!zo), then p.l7p reduces to the conservative shearing wave 
p.6p - (l3.8p . For 7^ 0, Gj^ — as t — /i — 00, which implies the eventual viscous decay 
of any shearing wave forced at a particular time fi. 

For the dynamo problem we assume that the velocity fiuctuations reach a stationary 
equilibrium after a finite time, long compared to tf and to an approximate viscous decay 
time, 1/ (fc^ji^). This formulation implicitly assumes nonzero viscosity, or else the random 
velocity variance would grow without limit and not equilibrate. 
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3.5. Kinetic Energy, Non-dimensionalization, and Homogeneity 
Define the volume-averaged kinetic energy as 



KE{t) = \ 



u 



x,y,z 



(3.20) 



where the angle brackets again indicate an average over the spatial coordinates. For this 
dynamo problem we adopt a dual normalization in the fluctuation forcing scale and in 
the resulting velocity scale, or equivalcntly the equilibrium kinetic energy: 



and E 



KE 



— — when 
2 



t~>t 



(3.21) 



Henceforth, all quantities are made non-dimensional by the implied length and velocity 
scales {i.e., forcing amplitude, time, magnetic field amplitude, viscosity, and resistivity). 
We further assume, for definiteness, that the expected value of kinetic energy (I3.2ip is 
equally partitioned between the horizontal and vertical velocity components in p.20p : 

1 



KE, 



£ 



KE\ 



(3.22) 



There are no cross-terms in KE because of the statistical independence of fz and 6^ in 
(|3.16l) . This partition thus gives separate normalization conditions for and Oz- We 
will see in Sec. Uthat both Fz and Oz must be nonzero for the shear dynamo to exist. 

For the solutions in (I3.17p . the kinetic energy density involves products of Fourier 
factors, with product phases ±0(^) ± i^(/z'), inside a double time-history integral over /n 
and /i'. The z average is trivially 1 for a barotropic flow with no z dependence in (f>. We 
assume the horizontal domain size L is large compared to the forcing scale, 1 /k±f. For the 
terms with summed phases, the x and/or y averages oi±2{kxfX+kyfy) are approximately 
if Lk±f :s> 1. (This could also be assured if Lkyf/2'K has an integer value as part of a 
discretization of the forcing; Sec. 15.21 ) Focusing on the remaining terms with differenced 
phases, we take an x average over phases ±{kx{t — ^) — kx{t — fi'j)x — ±Skyf{fi — ji'). 
After performing the z and y averages and substituting the forcing covariance functions 
p.l6p . the partitioned normalization conditions from p.22p are equivalent to 



Fz 



dA* / d^x' G,{t^ y)G,{t- ^i') 
Jo 

exp — 111 ~ ^'\/tf (exp iSkyf{ii~ii') 



= FzCz = 1 



d/i / d/i' G,(t-/i)G,(i-//') 



ki(<-/i) -ki (<-//') 



10 

exp 



klit-iJi) klit-ii') 



\li-fi'\/tf (^exp iSkyf{fi ~ fi')x ^ ^OzC 



= 1 



(3.23) 



which are independent of t as t — ^ oo. This defines the constants Cz and C± that then 
determine Fz and Oz - It will simplify the dynamo problem in Sec. l4.2l to renormalize the 
random forcing amplitudes by 



fl = Cl/^fz , 4 = C'foz , (3.24) 

whose corresponding expected variances are unity, F^ — CzFz ^ 1 and 0| = C±0± = 1, 
and the associated expected energies are KEz — F^/A and KE± — Ol/A. 

Cz and C_L are continuous, finite (if v > 0), and positive functions of S, v, L, tf, and 
the forcing wavenumber orientation angle Of, 

kxf = cos9f, kyf—smOf. (3.25) 
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Note that < 0f < 7t/2 is an up-shear tilt when S > 0, while tt/2 < 6/ < tt is down- 
shear. The extreme values 6/ — 0, tt {kyf — 0) are not of interest because there is no 
shear-tilting in (|3.18p and no dynamo in Sees. |3][ni 

We could proceed quite generally in all these parameters, but at the price of consid- 
erable complexity. Various degrees of simplification are available in different parameter 
limits, e.g., if the domain is large (as already partly assumed in Lk±f 3> 1), — >■ oo, 

5 0, or tf — > 0. The simplifications arise from being able to isolate and integrate 
over one or more of the factors in p.23p while approximating the time arguments of the 
other factors as fixed at the importantly contributing times insofar as they are varying 
relatively slowly. 

Among all these parameters, the simplifying limit that seems most physically general 
and germane is large L, with provisionally finite values for the other parameters. For the 
rest of this section and Sees. I4l l4.31 we follow this path, and in Sec. [5] some additional 
and alternative limits are discussed. On this path we isolate the spatial average factor 
in p.23p by the x-averaging operation explicit and integrating over its time argument, 

6 = fi — fj,', asymptotically over a large interval, while setting fi fj,' for the other factors 
(because the spatial average factor is small everywhere that S is not). Thus, 



dS ( exp 



iSkyfSx 



1 



OO 
— OO 



.L/2 




/ ds exp 


iSkyfSs 


/-L/2 





dS 



SkyfSL 



SkyfSL 



27r 



SkyfL 



(3.26) 



The final step on the second line is based on the asymptotic integral of the sine integral 
function. Si [mathworld.wolJram.com). To achieve this approximate isolation from the 
viscous and forcing-correlation factors, we assume LkyfS/v, LkyfStf 3> 1, along with 
the previous assumption for averaging, Lkyf 3> 1 This is not the distinguished limits of 
small S or tf in a finite domain (Sec. 15. 2p . although when taken successively following 
(|3.26|) such limits are well behaved (Sec. 15. ip . The relation p.26p can equivalently but 
more compactly be expressed as 



exp iSkyf{^ — ^')x 



(3.27) 



with Cl — 2TT/SLkyf. 

Inserting (pT27l) into ((3?23l) yields 



C 



Cz — <-^L^ _ _ 

Al = dpGlip) , Al = dpGUp)kljip) . (3.28) 

After the normalizations (|3.2ip - (|3.22p and the large L approximation yielding (I3.27p . 
the non-dimensional parameters of the BSD model are 5, v, tf, and 9f, plus other 
quantities related to B defined in Sec. HI There is no dependence on L. 



As an aside we examine the ensemble-mean local velocity variance, £ 



v?{x,y, z,t) 



which is different from the domain-averaged 28 KE . From p.l7p and the covariance 

properties of the random force p.l6p . e.g., the vertical velocity variance has the expected 
value at late time, 

/>oo 

dp I dp' G^{t — p)G^{t — p')Fz exp ~\p~ p'\/tf cos[Skyfx{p — p')]. 
Jq L J 

(3.29) 

This variance is independent of t because nonzero viscosity renders u stationary. It is 
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independent of y and z, i.e., homogeneous in these coordinates. But the local variance is 
not in general homogeneous in x. In the limit v ^ oo, the integrals can approximately be 
evaluated (as discussed more fully in Sees. 15.1] and 15)) to yield a constant value equal to 
= FzCz in (I3.23|) . For finite viscosity the peak variance is at x = 0, and it decreases 
with \x\ on a scale ~ l/{Skyftf); this can be seen by taking the limit of small tf where 



2Fbtf 



1 + (StfkyfX) 



AnGl{t-n). (3.30) 



Homogeneity is thus restored for small S or small t/, although these limits are formally 
incompatible with the approximation underlying p.27p . which is therefore to be under- 
stood as a horizontal average over a region that encompasses any variance inhomogeneity. 
The fundamental source of forced shearing-wave inhomogeneity is the special zero value 
of the mean flow SxSy at x = 0: the phase-tilting rate SkyfX increases with while 
the forcing correlation time tf does not depend on x. Homogeneity holds for v ^ oo 
because the forced shearing waves have non-trivial amplitude only for = (l)f, i.e., no 
phase tilting. 

An amelioration of the inhomogeneity magnitude results from the dynamical freedom 
to add a random forcing phase r(/i) to (|3.18p : e.g., a model for r is a 27r-periodic random 
walk with correlation time t,,. Inhomogeneity is eliminated if t^ — ^ 0, but it still occurs 
with finite tr. A broader posing of the BSD problem is for a family of mean fiows with the 
same mean shear, i.e., V = C/,e^ + (t4 -I- ^(a; — a;*))ej^, and a corresponding modification 
of the forced shearing- wave phase p. 181) to (j>{x,y,t; n) = k^it — ^){x — x^) + kyf{y — 
y^.) — k±fV^,{t — fi) + 0.5SU^,{t — fj,)'^+r{ij). An expanded-ensemble average over V, and 

over x» in particular, restores homogeneity in x of £ for general parameters, which 

thus is a corollary of translational and Galilean invariances. These generalizations in r 
and Y do not change the dynamo behavior in anything except the shearing-wave phase, 
which does not appear in KE or the BSD (Sec. 14.21 et seq.), so we now drop further 
consideration of them. 



4. Magnetic Induction 

Write the induction equation (I2.8P as 

VB X {ux B) + SB^ey + rjV^B . (4.1) 

To simplify matters, we note that the induction equation is linear in the magnetic field. 
Therefore, for a barotropic velocity field u{x,y), the electromotive force does not give 
rise to any mode coupling in z. We pose the dynamo problem as exponential growth 
of the horizontally-averaged (i.e., mean) horizontal magnetic field with an initial seed 
amplitude and a single z-wavenumber kz, 

(^BA_y'" ^Re{B{t)e''''^} . (4.2) 

Thus, both kz and the initial mean field, 6(0), are parameters of the problem; without loss 
of generality, we can take |0(O)| = 1 as the non-dimensional normalization of B. Because 
we are interested in dynamo behavior with exponential growth, this normalization choice 
does not affect the resulting growth rate. We then define Ob as its initial orientation 
angle: 

B4O)^cos0B, By{O)^sm0B. (4.3) 
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Because f{t) is a stochastic variable, the more precisely stated dynamo problem is expo- 
nential growth of mean magnetic energy looking over many realizations and/or 
long time intervals. 

Because there is no Fourier mode coupling in z, we can assume the entire magnetic 
field has only a single kz , and the application of the gradient operator is simplified to 

V = + ikzBz . (4.4) 

We only need to solve for the horizontal component of B, i.e., B±, and obtain 
diagnostically from the solenoidality condition, 

B.^-^. (4.) 

For the mean field (B±) , there is no associated vertical component. The horizontal 
induction equation from (14. ip is 

VBi_ = -{u • V) Si + (B_L • Vi) ui_ + SB^By + riV^Bi_ . (4.6) 

Because it is enough to focus on the horizontal components of B, we henceforth drop 
the subscript _L and interpret all vectors a as horizontal unless indicated otherwise by a 
subscript: a 3D vector will be {e.g., V3). 

The non-dimensional parameters in the ESD associated with the magnetic field are , 
r], and 9b', these are in addition to the dynamic parameters listed at the end of Sec. 13.51 

4.1. Magnetic Fluctuations 
Decompose the horizontal magnetic field into fluctuation and mean components, 

B{x,y,z,t) ^ 6B{x,y,z,t) + Re{l3{t)e'''''''} . (4.7) 

For consistency with (|4.2p . we specify that (SB) — 0. We evaluate the vertical com- 
panion field SBz by (|4.5p . Because (|4.6p is linear in B, we see that SB will have the same 
vertical phase factor as the mean field; i.e., we define its complex coefficient b by 

SB = Re{b{x,y,t)e'''^'} . (4.8) 

By assumption the ESD contains only a single phase component for the horizontal 
magnetic fluctuation field b(x, y, z, t) determined from the horizontal forcing wave number 
kf (through its shear-tilting phase <j) in p.lSI) ') and the vertical wavenumber k^ of the 
seed mean magnetic field. Its induction equation from (|4.6p is forced by the stochastic 
shearing waves and the horizontal mean magnetic field, i.e., 

VSB ^SF + SeySB^ + r/VlSB , (4.9) 

where the curl of the fluctuation electromotive force dF is 

SF{x,y,z,t) = -uA/^bY'"' + (^B^"'^ v) u 

= - UzRe {ikzBe'^'^] + (Re {Be'''''] • V) tt . (4.10) 

There is no contribution from — (m • V) (B) because (S) has no horizontal gradi- 
ent. One can view the ESD fluctuation induction equation (j4.9p for SB as a first-iteration 
approximation to the full MHD induction in the presence of u and \B) ; i.e., it is 
a projection of MHD onto a magnetic field with only the shearing-wave phase and a 
horizontally uniform component. 
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This simplified equation for h is the heart of the quasi- linear ESD theory {i.e., linear for 
magnetic fluctuations, nonlinear for the horizontal mean). The quasi-linear simplification 
can be rigorously justified only if |6| <C in which case all higher harmonics of the 
phases in b will be negligibly small compared to the primary phase; this condition is 
met in the limit 77 — 00, i.e., vanishing magnetic Reynolds number (Sec. [5]). In the 
next section we will see how the spatially-averaged induction from the shearing waves 
induces dynamo growth in 3. This quasi-linear theory is formally incomplete when the 
preceding justification condition is not always well satisfied by its solutions. Nevertheless, 
they correspond to th e shear dynamo foun d i n 2+D and 3D simulation s for a fairly broad 
range of parameters ( Yousef et al. 200830; iHeinemann et al. 20116[ ). so we infer that 
the ESD provides a cogent explanation of the dynamo process even beyond its rigorously 
derivable limit. When u variance is inhomogeneous fSec. 13. 5|) . SB variance will be so as 
well. 

Using the shearing wave solution (I3.17P and the mean field expression in (j4.7p and an 
analogous vertical phase factor decomposition for ^i^ as for SB in (|4.8p . we evaluate the 
fluctuation forcing term as 



Fb{x,y,t)= dfiG,{t-fi)[-ik,B{t)Re{f,{fi)e"^^^^} 

(fc(i-M).B(i))Re{6,(A.) e'^(^)} 



X k{t — jl) 



(4.11) 



Pro tem we do not yet use the renormalized forcings p.24p but will do so in the next 
section. The two right-side lines here are, respectively, from the two terms in the second 
line of (|4.10p . The magnetic fluctuation Fourier phases are thus ± </>(/i) -f- k^z where is 
the shearing wave phase in (|3.18p . 

We can write the solution of (|4.9p for b analytically. Again utilizing the vertical phase 
factorization (|4.8p . we have 

t I'X 

dA / d/j, Grf{t — fi, X — fi)G^{X — fi) 
'0 Jo 

- ifc,S(< - A) • B(A) Re { M^l) e'^(^) } 
( fc(A - pt) ■ B{X) ) Re { d,{fi) e'-^^^) } ] . (4.12) 



b{x,y,z,t) 
Bz X k{t — fl) 



k^X-fi) 

Here we define the second-order, real tensor S by 

S{t) = I + StByB^ , (4.13) 
with / the identify tensor, and the resistive decay factor (another Green's function) by 

Gri{t, X, fi) — exp — rj / dpk^{p — ^) 



exp 



1 I ACkiiO 



(4.14) 



with fc|(f) = fc^(i) + k^. Again, for compactness we write this as G,,(i — /i, A — /i). Thus, 
in the quasi-linear ESD, SB is an induced magnetic shearing wave arising from u and 

^> • 

4.2. Mean Field Equation 
The governing equation is the horizontal average of ([4.6^ : 

x.y j \ x,y j \ x^y 



dd B 



= \ -f B 



ey-vKi B 



(4.15) 
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where 



Fb y^\z, t) = (^-{u- V)b' - {u,d,)b' + [b' ■ V)w ) 



(4.16) 



Because of the horizontal average in the ESD mean-field equation, there is no repre- 
sentation of any spatial structure associated with wave-averaged inhomogeneity in the 
electromotive force curl (Sees. [331 and l4.1( ). 

The induction forcing itself depends linearly on B through b in ()4.12|) . where it enters 
in a time- history integral. So ()4.15|) is a linear integral-differential equation for B{t), for 
which no general analytic solution is known. Instead, we evaluate the expression for Fjg 
below and obtain a double-time integral, second-order tensor operator on 13(t) that we 
will solve numerically in general (Sec. l4.3() and analytically in certain limits (Sec. [5]). This 
yields a closed-form equation for the mean magnetic field amplitude as a function only 
of the forcing time histories, fz{t) and dz(t), and the parameters kf, S, rj, and v. 

As with the b solution in the preceding section, the derivation for ( i^g ) is rather 

elaborate. It involves substituting the shearing wave solution (I3.17P and the magnetic 
fluctuation (|4.12p into (j4.16p and performing the horizontal average by identifying the 
zero horizontal phase components and applying p.27|) : these details are in Appendix lAl 
If we again define a vertical Fourier coefficient F13 , as in (|4.2|) , the result is 



Notice that the forcing helicity -ff (/i) from p.lSp plays a prominent role. With the so- 
lutions in Sees. [Sjni we find there is only transient algebraic growth in B{t) {i.e., no 
dynamo) when the forcing helicity is zero. Therefore, there is no dynamo if either fz 
or 6z is zero. In fact, the induced magnetic fluctuations from a horizontal velocity held, 
forced by dz only, have no effect at all on B. 

Now simplify Fq and the B equation by the forcing renormalization p.24p augmented 
by the following related quantities: 



-^^ = ^1/]^ ^^ = ^Re{/t*(^)6t(^)} = ^i?, Gl = j^G.. (4.18) 



FB{z,t) 




(4.17) 



-z 



With these the mean electromotive force curl becomes 




"'0 



I T^fi) kl S{t - A) • B{X) + ikz n^'in) Bz x k{t - ^) 



(4.19) 



An identity used for the flnal term is k{t — fi) ■ S{t — A) • B{X) — k{X — ^i) ■ B{X). 

After factoring the structure Re { • e'*^^^} from (|4.15p . the equation for the complex 
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amplitude B{t) becomes 
dtB = SB^e. -'-2 



t n\ 

dA / 

^0 



d^i Griit - M, A - ^l)Gl{\ - ii)Gl{t - m) { klS{t - A) • B{\) + 

ik,n\^l) X k{t - ^l) k-^{x- ti) + k-^{t- n) (fc(A- /x) • b(a)) } . (4.20) 

A final compaction step is to factor out tlie resistivity effect associated witlr tlie vertical 
wavenumber by defining 

B{t) = 66-"'='* . (4.21) 

This modifies (IT^ to 

dtB = SB^ey 

* ~ r 

dA / dfiG,j{t-i^,X-^)Gl{X-fi)Gl{t-fi)\P{iu)klS{t-X)-B{X) + 
Jo ^ 

ik^n\li) X k{t - m) ( fc"^(A - m) + '^"^(i - m) ) ( fe(A - A*) • ^(A) ) } , (4.22) 

where G,, is the resistive decay associated with the horizontal wavevector, defined analo- 
gously to G^ with fc3((^) replaced by k{(^) in (|4.14p . i.e., factoring out the decay associated 
with kz, 



(t , A, A* ) = exp - r]kl {t ~ X) G^ (i , A, /x) 
The functional form of (|4.22p is 



dtB^L- Bit) 



dA J{t,X)-B{X), 



(4.23) 



(4.24) 



where L and J are second-order tensors. This ESD form differs from the common ansatz 
Jl.ll) by the time-history integral, but it does fit within the formal framework analyzed by 
ISridhar fc SubramanianI ( 2009.) for velocity fields whose dynamical origin was unspecified 
(in contrast to our particular case of shearing wave velocities). We show in Sec.[5]that the 
common ansatz is recovered in our ESD theory in the limit of 77, — >■ 00. The definitions 
of the L and J tensors are 

Jmn{tA)=- I d^lG„{t~ ^1,X- ^i)Gi{X~ fl)Gl{t- fl) 

Ja 

THp)kl {s^mn{t-x) + ik,nHf^) 

( k-^{X + k-\t - ^i) ) kt{t - n) kn{t - A) ezim ] (4.25) 

for horizontal indices, {m,n,£} = {x,y}, and the usual Kronecker delta and Levi-Civita 
epsilon tensors. L contains the background shear effect on B, while J contains the mean 
electromotive force resulting from the random barotropic forces and induced magnetic 
fluctuations. Smn = Smn + S {t — X)6mySnx IS as defined in (I4.13|) . 

The ESD (j4.20p is invariant with respect to several sign symmetries in the forcing, 
wavenumber, initial conditions, and mean shear. Because the random forcing amplitudes, 
fzit) and 6z, are statistically symmetric in sign, a change of sign in either one implies 
O — 'H^, J-"^ -f-)- J^'' , and the statistical distribution of B{t) will be unchanged. In 
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addition there are the fohowing invariances for particular reahzations of the ESD: (i) 
{h, H^) o -{k,, Ht); (ii) kf o -fc/; (iii) - B; and (iv) {S, V.\ k,f, B^) o 

- (5, n\ k,f, B,) with (kyf, By) ^ {kyf, By). 

Because the ESD is a quasi-hnear theory based on Fourier orthogonahty in kz and 
kf, it satisfies a superposition principle; the full MHD equations (|2.ip - (|2.4p do not allow 
superposition, of course. The functional form of the superposition is a generalization of 
gill) and KM^ : 



(z, = ^ Re i) exp —rjklt + ikzZ | , 

dtB{kz,t) = L{kz)-B{kz,t) + dX |^|]j(fc,,fc/,t,A) j 



J(fc.,A). (4.26) 



The random force f{kf,t) in p.9p is assumed to be statistically independent for each 
k f component with whatever normalization is chosen in place of the single-component 
normalization (|3.21l) . 



4.3. Dynamo Behavior 

A numerical code has been written to solve the ESD in (|4.22p . Its algorithm is described 
in Appendix |B] As expected from the 3D and 2+D full PDE solutions, a dynamo often 
occurs when S and 'H(t) are nonzero. We now demonstrate a typical dynamo solution, 
deferring the more general examination of the ESD parameter dependences until Sec. [6l 
after first obtaining analytic solutions in Sec. [5] in certain limiting cases. 

An illustration of a random realization of the forcings, velocity variances, and helic- 
ity time series is in Figs. [T][2] These are for a case with moderately up-shear forcing 
wavenumber orientation {9f = 7r/4), moderately small correlation time tf = 0.1 and 
viscosity v = 0.1, and intermediate mean shear rate {S — 1). The amplitude normal- 
izations from p.22|) are evident, as is the vanishing of the time-averaged helicity. Be- 
cause tfV <^ 1, the time scale of the velocity fluctuations is controlled primarily by 
the viscous decay time modified by the shear tilting in the kx{t): in p.l9|) the initial 
exponential linear decay rate, ly = 0.1, is at first slowed as kx passes through zero at 
t = l/StanOf = 1 and then augmented toward a exponential cubic decay with a rate 
coefficient w {vS^kl^/S)^/^ = 0.26. 

To obtain a dynamo in (|4.22p . the vertical wavenumber kz must be small but finite; 
we show below that this is true for general parameters. With kz — 0.125 and moderately 
small r] — 0.1, the time series of the mean magnetic field component variances are shown 
in Fig.|3]for the same realization of the forcing and velocity as in Figs.[T][2] There is evident 
exponential growth in both components of B{t), i.e., this is a dynamo. If we make an 
exponential fit over a long time interval with \B\ cx e'^t, we obtain the same value of 
7 ~ 0.03 for each component. B{t) also manifests a stochastic variability inherited from 
the random forcing, and its fluctuations about the exponential growth exhibit power even 
at much lower frequencies than are evident in the forcing and velocity time series. 

The magnitude of By is larger than of Bx in Fig. [S] This is a common behavior for 
magnetic fields in shear flow. A partial and somewhat simplistic explanation is as a 
consequence of the first right-side shear term in (|4.22p . A simplified (non-dynamo) system 
with arbitrary forcing R{t), 



dtB ^ SBxey + r{t) , B(0) = Bq, 



(4.27) 
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has the solution, 

B{t) = Bo + dt' R{t') + Sey(^BMt + J^^ dt' dt" r^{t")^ . (4.28) 

The last term oc Sey will make \By\ ^ \Bx at late time for most R{t). This anisotropy 
effect carries over to the ESD but also involves further right-side B coupling absent in 
(j4.28p : a coupled explanation for the anisotropy in dynamo solutions is made in Sec. 
15.11 The initial condition Bq is usually not dominant in (I4.28P at late time. The initial 
condition is even less important for B{t) in Fig. |31 which is obtained with 9b = 7r/4; in 
particular, 6b does not determine the dynamo growth rate 7. 

b{t) (not shown) also shows exponential growth in its amplitude, with \by\ typically 
much larger than \bx\ for the same reason as just explained. b(t) has comparable time 
dependence to Uz{t) and u{t), as well as an additional resistive decay influence from 77 
and modulations by the exponential growth and slow variation in B{t). 



5. Dynamo Analysis in Limiting Cases 

5.1. L — > 00; T], v ^ 00 

The ESD in Sees. [5]|3]is based on an assumption that the horizontal domain size is large, 
L — ^ 00 (n.&., the average of a Fourier exponential in p.27p ). As a means of obtaining 
a more readily analyzed form of the ESD (I4.22p . we take the additional limit of 77 — >■ cx). 
This limit does not change the forcing amplitude nor the velocity field (Sec. which 
are independent of rj, but it allows an elimination of one of the time integrals in the 
expression for b in (|4.12l) and in the equation (I4.22p for B. It also makes the quasi-linear 
approximation rigorously accurate because it yields \SB\ ^ \B\ (as explained after (|5.3p ). 

The essence of the ry ^ 00 approximation is that first the order of integration in (|4.22p 
is reversed, 

dA / d^= / d^ / dA, 

Jo Jo Jt-n 

and then the A integral is performed by assuming that G,, is more rapidly varying in A 
than any of the other integrand factors and furthermore is nonzero only when A — )• t, 
i.e., t — X = 0{ri~^). We evaluate this approximation as 

d\G^{t~^l,\~^l)^-^^^^—^ (5.1) 

for all /i 7^ i (the integral is zero for fJ- = t) and set the A arguments of other factors in 
the integrand to t. With this approximation, the (i, 77)-limiting form of (I4.22p becomes 



dtB = SBxit)ey - - / d/i 
V Jo 



1 f ^ qlit-jj^ 
k^it-^i) 



klFHp.) Bit) + 2ikMHp) ''\2(f_^ ^)^^ '^(i - m) ■ m ] ■ (5.2) 

This is a purely differential equation for B{t)] i.e., it matches the common ansatz form 
in (jl.ip . viz., 

dtB= L- B{t) , (5.3) 

for the identifiable single-time, second-order tensor L{t) that contains a time-history 
integral in jj, over the random forcing. 
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An analogous simplification of the expression for b in ()4.12|) can be made, with the 
result that b cx This gives the important analytic result that the quasi-linear ap- 
proximation to (j2.3l) is asymptotically convergent as 77 — > oo; the higher harmonics of the 
shearing-wave Fourier phase (±m(/>, m > 1) generated in b by the fluctuation electromo- 
tive term are 0{ri~"^), hence negligible compared to the mean-field term proportional to 

(^bY'^ in Kim . 

Numerical solutions of (j5.2p exhibit dynamo behavior similar to the example in Sec. 
14. 3[ and the parameter dependences for 7 are similar to those described in Sec. [6] for the 
general BSD. In particular, 7 is small here because 77 is large, in contrast to the "fast 
dynamo" limit where 7 becomes independent of rj {cf., Fig. I10|) . 

To obtain further analytic simplicity we can take a sequential limit of (j5.2p as — ^ 00. 
As with the ry limit, this selects an integration time /i ~ i, where the viscous decay factor 
is integrated out by the approximate relation for large t, 

\^,GUt~^l)^ 1- or f d^^Gl^it^^i)^!, (5.4) 

2;^ Jo 

utilizing the renormalization relations in (|3.28p and (|4.18p . The {L, 77, z^)-limit mean-field 
equation from (j5.2l) is 

dtB = SB^{t)ey ~ i klJ'Ht) B{t) + 2ik^'H\t) (e, x kf) kj ■ B{t) , (5.5) 

after using fc^(O) = fcj = 1 from p.21|) . In the tensor representation ()5.3|) . L{t) is defined 
for (I53|) by 

0\ klJ'Ht) fl 0\ 2ik^n\t) /cos 6>/ sing/ sin^^/ \ 
^l^lOy rj 1^0 ly rj \ — cos^ 9f — cos 9 f sin 9 f J ' ^ ' ' 

after a substitution for kf from p.25p . All of the forcing time history in the coefficient 
tensor L{t) has now disappeared. The history integral also disappears in the companion 
b formula derived from (j4.12l) . Furthermore, there is no remaining dependence on v in 
(|5.5p because J-"^ and are 0(1) quantities by the KE normalization in p.2ip and the 
forcing renormalization in p.24p and ()4.18p . Large 77 and v values lead to momentum 
and induction equation balances with negligible time tendency terms and negligible shear 
tilting in k{t) because (j> (pf and k(t) — >■ k/. 

We now consider two further limits in the forcing correlation time t / that yield analytic 
expressions for 7. 



5.1.1. Steady Forcing 

Suppose the forcing values taken from the random distributions in Sec. 13.31 but are 
held steady in time; this is a limit based on the physical approximation that the forcing 
amplitudes change more slowly than the inverse growth rate for the dynamo, jtf ^ 1. 
In this limit ()5.5p - (|5.6p has its L independent of time, hence there are eigensolutions with 

B oc e^K (5.7) 

The eigenvalues of L are 

r = - Mj-t ± / 2ifc. sin^ 9fH^S 
V \ V 
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The dynamo growth rate for total mean field 13 is defined as the largest real part of V 
plus a correction of — rjk^ from the transformation in (|4.21|) : 

The first term is negative and the second positive. A dynamo occurs with 7 > if there 
are both forcing helicity and shear and if is small enough but nonzero. With 5 = 0, 
there is no dynamo. For \S\ above a critical-shear threshold value, 

^- = ^ — w\ — 

7 increases with S*, asymptotically as when the other parameters are held constant, 
and 7 decreases with 77 as I/77. For given 5, there is a lower threshold value for 77 to 
have a dynamo. Nonzero forcing helicity is necessary for a dynamo, but its sign does not 
matter. 7 = for kz — Q, and 7 < for fc^ large. Within an intermediate range where 
7 > 0, the optimal and its associated growth rate are 



opt — 



|gHt|sin^g/ \^'\ (\Sn^\sin^ Of^ 



( 27|5Htpsin^[g/] ^ /27^Htpsm4[^^ 



Use r?2(r, + 77-1 J-t)j ~ (, 256r;3 j ' ^^'^^^ 

where the approximations are based on neglecting F'^)/irj^ . The optimal kz decreases with 
increasing rj. (In a general MHD simulation with fixed (5*, 77, v) values, all kz are available, 
and the ones supporting a dynamo will emerge in the evolution.) The vertical forcing 
variance reduces the dynamo, while the forcing helicity amplitude \'H'^\ enhances it. 
J^t enters (|5.6p and (|5.9p exactly as an enhanced resistivity; however, the effect is small 
as 0(77^^) when J^t = 0(1) in this large 77 limit. This is an anisotropic turbulent eddy 
resistivity acting on the mean field in the direction perpendicular t o the shear plane as a 



result of the shearing-wave vertical velocity (Parker 1971; Moffatt 19781) . The horizontal 
force / acting by itself has no effect; it makes F — =0, hence 7 < (no dynamo). 
7 is largest where kyf is largest at Of = 7r/2; in Sec. |6]we show that 7 is usually larger 
for Of < 7r/2 (Fig, ni]) because of a dynamo enhancement by the shear-tilting Orr effect 
when V < 00. k^f does not explicitly enter the formula for 7 in the present case. 

The system (I5.5p - (l5.6p in its steady- helicity l imit is a close a nalog of the so-called 



ihe system (l5.op - (lo.b|l m its steady- nelicity l imit is a close a nalog ot the so-called 
alpha-omega dynamo for galactic disks IParkei ( 1971 ) ; Kulsrud ( 201Clf ) . Using a mixed 



x.y 

notation from these two sources and assuming a vertical structure { B ) cx e^'^'^ , an 



= ( ) ■ (5.12) 



ODE system analogous to (j5.3p results, with 

' 77fc: 

Vl —rjk 

For constant a and fi, its eigenvalues are 

T"" ^ ^ klJ^ ± iikzan)^/\ (5.13) 

The correspondence with (j5.8l) is evident with appropriate identifications between (a, rf) 
and {'q~^'H\ S, 77 + rj^^J-^). However, the ODE systems are not isomorphic except in 
the special case of kxf = in (|5.3|) . Thus, in the steady- forcing ESD, the shear S plays 
the role of fi, helical forcing plays the role of a, and J^^ plays the role of a turbulent 
eddy resistivity that augments the effect of 77. 
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The physical paradigm in this paper is random forcing. Therefore, even if the forcing 
is steady in time, it is taken from a random distribution, and we can ask what the 
expected value is for B {i.e., having factored out the resistive decay in (|4.21l) . which is 
not dominant for small kz). To answer this we now neglect the turbulent resistivity by 
J^'' , which is shown above to be a small effect for large rj. The eigenvalue (|5.8I) of the 
tensor (j5.6l) is for a particular forcing value, which we now generalize to an ensemble 
distribution, 

r(e) =±7(l + is), j{e) = ^ESsinOf > 0, (5.14) 
with a composite parameter that is a rescaled helicity forcing. 



e = 



Sr/ 



= E^s. 



(5.15) 



£'2 is the magnitude of e, and s = ±1 is its sign. Consistent with the Ornstein-Uhlenbeck 
process for the forcing amplitudes fSec. [?75)) . e has a Gaussian probability distribution 
function. 



V{e) 



1 



exp 



Vde^l, 



(5.16) 



with an expected variance Eq. Utilizing £ = 0.5 Fz] 0\ = 0.5 from the remark after 

(|3.24p . we obtain 

2kl 1 



el 



where the arrow indicates substitution of from (|5.1ip . We analyze the dynamo so- 
lutions with general Eq, but for large r?, Eq is expected to be small. After a large elapsed 
time ig) the dynamo solution is dominated by its leading eigenmode with Re { F} = 7 > 
for any E 0. Neglecting the decaying mode, we write the late-time solution in vector 
form as 



I Bx(e,te)\ ^ -yt / r , 1 . • r , IN /(I + islT/S* + iecos6'f sin 
Co is a complex constant determined from the initial condition. 



Co 



1 



\/2E{i + is) 



(1 + is)7/S' — iecos6'/ sin6'/ 
1 — ie cos^ Of 



(5.18) 



(5.19) 



With (|5.16|) and ()5.18p . we can evaluate the expected value of any property of 3{te) and 
its corresponding distribution D with e; e.g., for the mean-field vector magnitude, 



£ 



1^1 {te) 



\B\{e,te)V{e)Ae 



D[\3\]<ie. 



(5.20) 



Figure m (left panel) shows the distributions D for the vector magnitude and for the di- 
rectional component magnitudes for a small value of Eq- These distributions are smooth, 
positive, symmetric in s, and peak at intermediate e/eo- B™'^ and the component mag- 
nitudes are growing exponentially with time. We can fit this with a cumulative growth 
rate, Y'^^ = logi?'"™", which we know from (|5.14|) will scale as Sy^sinOf. For this 
value of eo = 0.1, is smaller than \By\, with a ensemble-mean ratio of 0.78. For the 
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leading eigenfunction in (|5.18p . the anisotropy ratio is 

\Bx\ (1 + is)Esa\9f cos 9f sin 9f 



%/2(l-i£cos2 9f) 



(5.21) 



For small E, the ratio tends to E sm9f / \/2, which is small; this is consistent with the 
anisotropy in Fig. [31 For large E, the ratio tends to | tan^j^j, which can have any value. 
What is the ensemble-mean magnetic field? Its magnitude is 



B{s,Qrie)\de 



D[B]de 



(5.22) 



Again this is evaluated with ()5.18|) . We find that it too exhibits exponential growth, so 
we fit a cumulative growth rate, 7'"^°"(te) = t~'^ logB™^"" > 0. But the ensemble mean 
growth is smaller than the ensemble r.m.s. growth, i.e., Y"''^"''^ < 7''™*. The reason is 
illustrated in Fig. |4] (right panel) for the distributions of two components, £'[Re | /B^; | ] 

and _D[Im | | ]. Their amplitude is comparable to the magnitude distributions in the 

left panel, but they are oscillatory in e as a result of cos[7te] and sin[7te] terms in (|5.18p . 
So the expected value from integration over e is small, although not zero. For Fig. |4l 
^mean ^ o.073B™'^ and 7""=''" = 0.767™'^ 

These relations are not sensitive to the initial condition S(0), although it does influ- 
ence the partition among the real and imaginary parts of 13{te). There are the expected 
dependences of larger 7 with larger S and Eq and with 9f closer to tt/2, as in ()5.9p . With 
larger te the expected values are dominated by the farther tails of the D distributions, 
with slowly increasing 7''™*(te) and 7™'^""(te) associated with larger 7(e) in the tails 
(Fig. [5]). Even though larger e values are less probable in P{e) in (j5.16p . they do have a 
more than compensating stronger dynamo growth rate that emerges after long enough 
time. Because the discrepancy between 7'''"* and 7™^'^" persists even in the P(e) tail, the 
ratio i3™eora^j^rms ^gcreases with te exponentially. The steady- forcing dynamo does not 
become independent of tg &s te ^ 00, in contrast to the finite-^/ dynamo, in particular 
the small-</ dynamo analyzed in Sec. 15.1.^ 

5.1.2. Rapidly Varying Forcing 

The limiting forms for the ESD equation, (|5.2I) and (|5.5p . are also analyzable in the 
opposite limit of — » by means of a cumulant expansion of a linear, stochastic, ODE 
system (jvan KampenI 120071 Chap. XVI). For a stochastic vector A{t) governed by 



dtA={Lo + U{t)) -A, (5.23) 
with the tensors Lq independent of time and Li{t) a random stationary process with zero 



expected mean and finite variance, the expected value £ 
deterministic ODE system. 



satisfies the approximate 



dfS 









A 




/ ^ 






Jo 



L,{t) U{t-t') 



dt' 



E 



(5.24) 



with the dots indicating neglected higher-order cumulant terms. The system (|5.24p has 
a time-independent matrix; hence, it has eigenmodes with exponential time dependence 
with growth rates given by the matrix eigenvalues. The solution formula for 8 A (t) is 

called a time-ordered exponential matrix, and it has a non-terminating series expansion 
with the leading terms as indicated here. The basis for the approximate neglect of the 
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higher order terms can be taken as the vanishing of £ Li^{t) Li{t — t') except as |i — 1'| — >■ 
0. In the present situation with large i/, this is equivalent to short correlation times i/ — >■ 
for the random forces, fz{t) and Oz(t), with Stf and S/u <C 1 to be able to neglect 
higher-order products of Lq and Li in deriving (j5.24p . 

We apply to with A = B exp k^T^t/T] with the following tensors: 



S 



^ _2ik^nHt)_ /cos 61/ sin 6*/ 
^ rj \ — cos^ f 



cos 



/ 

/ sin Of 



(5.25) 



This is a second-order, complex system. We have made one ad hoc simplification here, 



= 0.5 from (|4.18p . and then 



viz., replacing T'^{t) by its expected value, J"q = f J^t 

factoring its decay effect on B analogously to (|4.2ip . The motivation is to simplify the 
analysis. We already understand as an eddy resistive damping. This role is played with 
qualitative fidelity by retaining only its mean value, and anyway for large rj it is only 
a small increment to the ordinary resistivity. The result for (j5.24l) is very simple with 
(|5.25p because = independent of its time-variable prefactor, and the eigenvalues of 
Lq are zero. Hence, again after restoring the resistive decay factors, the growth rate for 



B 



IS 



7 = 



2^1 '^^^'-^ 



(5.26) 



i.e., in this (77 — >■ c», </ — > 0) limit there is only resistive decay of the expected value of 
the mean magnetic field, weakly augmented by the eddy resistive effect. 

We cou ld continue the cum ulant expansion for B and (|5.23p to higher orders in Stf 
and S/v (jvan KampenI 120071 ) . seeking growth in the ensemble-mean, large-scale field, 



B 



, but its 7 would be small in these parameters compared to the growth in the mean 



magnetic variance, £ 



\B\ 



To obtain a dynamo result for the latter, we instead apply 



(|5.23p - (|5.24p to the fourth-order real covariance system derived from (|5.2p for the vector, 
A - ( \By\^, Re[B:By], ImiB^By] ) x exp[2kl4t/v] , (5.27) 



again factoring out the mean eddy resistive effect with the simplification J^^{t) 
0.5. The associated tensors are defined by 



(0 
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1 











^0 








0/ 







sin^6'/ 2cos( 







'/ sin 6'/ 



2sin^ 6f 
2cos2 Of 
-2 cos Of sin( 




.(5.28) 



The expectation value in (|5.24p applied to Li{t)Li{t — t') acts entirely on its scalar 
prefactor in (|5.28p because its matrix factor L\ is deterministic and time-independent. 
We evaluate the corresponding scalar prefactor that arises in (|5.24p as 



4^ 

jy2 



dt' . 
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Tracing backwards through the forcing relations (|3.15|) . p.24p . and ()4.18p . we derive 



E 



= 0.5£ 



t|2 



E 



^t|2 



exp 



(5.29) 



utihzing the fact that the real and imaginary parts of and 62 are independent, station- 
ary processes each with an exponential correlation time tf as in p. 161) . After performing 
the time integration with this expression, the value of the preceding prefactor is 



2kltf 



n 



t2 



k^tf 



(5.30) 



because E 



0.5 F2 ] 0\ = 0.5 from (I3.24p . This completes the specification of the 



deterministic, time-independent matrix in (|5.24l) for the covariance system as 



L = SLl 



(5.31) 



We evaluate its eigenvalues F analytically from det[L — Tl] = 0, which is a fourth-order 
polynomial equation. We can factor a F = root, leaving a third-order system with the 
reduced form of F^ + pT = q for coefficients p (x S and q oa S^. With a simplification 
provided by the prefactor (j5.30p being small compared to S, we can neglect the p term 
and obtain the approximate solution. 



-,1/3 



2fc^5'^ sin^ 0f tf 



1/3 



(5.32) 



This approximation is consistent with finite S, small tf and k^, and large ry; recall that 
we also assume Stf, S/v <C 1 for the leading order cumulant approximation (j5.24p . 
The three solutions (j5.32p are one with real, positive F {i.e., a dynamo) and a complex 
conjugate pair with Re[F] < 0. We divide the positive eigenvalue F by 2 and restore the 

resistive decay factors to obtain the growth rate for the r.m.s. value of the mean field, 

1/2 



{£[\B?]) 



klS^ sin* 



^/ tf 



477^ 



1/3 



(5.33) 



A dynamo can occur with 7 > if there are both forcing helicity and shear and if k^ is 
small but nonzero; this behavior is the same as in the steady-forcing dynamo (|5.9p for 
this same limiting ESD system ()5.5|) . as well as for the general dynamo in Sec. El In this 
limit of small correlation time with zero mean helicity and finite helicity variance, the 
expected value for the mean field B does not grow, but the expected value for the mean 
magnetic energy 3^ does. The steady-forcing dynamo also has a much smaller ensemble 
mean than r.m.s. fSec. 15. 1^1]) . 

Besides the leading eigenvalue (I5.32[) . we can obtain the associated eigenfunction for 
the matrix (|OT|) . With the same approximation of a small prefactor for /.^^ , we derive 
the following for the expected ratio of component variances. 



2 sin'' e 



E 



F 

2fc2sin'*6lfi 



IS,, 



S'772 



2/3 



\B„ 



(5.34) 



Thus, the streamwise mean magnetic energy is small compared to the transverse energy 
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in the present hmit with transient forcing, smaU kz and tj, and large rj. The small 
ratio is also consistent with the previous example of dynamo behavior with more general 
parameters in Fig. |31 as well as with the steady- forcing dynamo in Sec. IS.l.T] when Eq is 
small. 

As with the steady forcing (jS.lip we can optimize the growth rate in kz- 




(5.35) 

The parameter tendencies here all have the same signs as with steady-forcing and with the 
general ESD (Sec. [5]), but the exponents are different in the two tf limits. In particular, 
the optimal growth rate dependences are 

7-5" \\n\\r]'^/^klft]/^ astf^O 

7-- ^2/3||Ht||2/3,y-ifc4/3^o asi/^oo, (5.36) 

where the norm symbol || • || denotes the r.m.s. or mean magnitude as appropriate, and 
we have formally restored the helicity variance factor WHW for emphasis. In both cases 
the growth rate 7 is vanishingly small as ry — )• 00, S* — 0, WHW — > 0, or Of — )■ 0,7r, and 
for the short correlation time case, 7 is small as — >■ 0. For non-limiting values of the 
parameters, however, 7 is not small (Sec.[6|). We reiterate that there is no dependence of 
7 on 1/ in the limit v ^ 00, independent of the value oi tf. 

As with the steady forcing limit, an analogy exists between the fluct uating helicity ESD 

in (15.511 and a low-ord er ODE fluctuating alpha-omega dynamo ansatz (jVishniac fc Brandenbure 
ISilant evll2000[ ) (also called the incoherent alpha-shear dynamo). Therefore, from 



a historical perspective of astrophysical dynamo theory, we see that the ESD in (|4.20p 
provides both a theoretical justification for the alpha-omega ansatz, with an explicit 
characterization of the relevant shearing-wave velocity fiuctuations, and a generalization 
to finite Reynolds numbers (i.e., ri,i> < 00). 

In summary, these two different t / limits with analytic dynamo solutions for the large- 
(j], v) ESD (j5.5p show qualitatively similar but functionally different parameter tendencies 
in S, T], kz, and Of] anisotropy with \By \ usually larger than \Bx\', and an ensemble- mean 



magnetic energy, E 



much larger than the energy of the ensemble-mean field. 



IB 



These characteristics carry over to the more general ESD solutions in Sec. [6l 



5.2. Other Limit Pathways 

I \ x,y 

The preceding ESD derivation of (j4.22|) assumes kyfL ^ 1 to assure ( ea;p[i(0 -I- (j)')] ) ss 

and kyfLS min[t/, l/i^] S> 1 to assure (exp[i{(j) + <j)')] ) 7^ for selected time argu- 
ments of the phases 0(/i) and 0(/x'). The latter assumption yields p.27p . which is useful in 
simplifying the normalization condition (I3.23P for KE and compacting the ESD equation 
(|4.22p for B by reducing the number of time history integrals in the mean electromotive 
force curl (Appendix [X]). We prefer the physical rationale of this pathway based only on 
a primary assumption of large L, consistent with uniform mean shear and no boundary 
conditions, because it does not constrain the values of the other parameters that are 



The Elemental Shear Dynamo 



23 



physically more meaningful than L. The result is independent of L itself. The further 
ESD simplifications in Sec. 15.11 follow from ry, — >■ 00. 



However, this is not a unique pathway for deriving ESD equations that are essentially 
similar. In particular, neither of the limits 5 — >■ nor t/ — )• is problematic even though 
they appear inconsistent with the second assumption above. As previously explained, 
we do require v > Q for statistical equilibration of velocity fluctuations and kyf 7^ for 
nontrivial shear tilting and dynamo behavior. 



Shear tilting makes 4){t) in (13.71) or p.lSp a continuous function of time. When 5 = 0, 
(j) — (pf, and the average of the differenced-phase factor is ( exp[i{(f) — (f)')] ) = 1 for all 

time arguments. When 5 — > as a primary assumption, this relation is approximately 
true. We still require the weaker assumption about large domain size, kyfL ^ 1, to 

be able to neglect the summed-phase factors, (exp[i{cj) + 0')] ) . Even with these phase 

averaging relations resolved, further assumptions are needed to compact the electromotive 
forcing, and large rj and/or i> suffice. The outcome is equivalent to (j5.5p with dynamo 
solutions when 5' > 0. If instead the primary assumption is — >■ in combination with 
kyfL ^ 1, then the requirement on the average of the differenced-phase factor in the 
KE normalization is resolved with an approximate integral over the forcing correlation 
factor, exp[— l/i — ^'\/tf], in p.23|) . but this assumption is not enough to compact the 
electromotive force curl. Again this can be accomplished with additional assumptions of 
large t] and/or z/, leading to the equivalents of (|5.2[) with shear tilting and (jS.Sp without 
it. In neither of these limits is there a compact equivalent to the general ESD (I4.22p with 
finite 1] and v. Also, because the dynamo solutions of (jS.Sp have 7 small with S and t/, 
this derivation pathway is not as physically germane as the primary one in Sec. 15.11 



Yet another derivation pathway assumes finite L and spatially periodic boundary 
conditions in shearing coordinates with discretized shearing-frame wavenumbers with 
Ak — 2Tr/L. If the forcing is at one of the discretized wavenumbers at least in kyf, then 
the spatial average of the summed-phase factor vanishes. To accommodate continuous 
shear tilting in the finite Fourier series representation, the forcing amplitude time series 
is viewed as impulses at discrete times, im = + mAt, At = 2tt / Sky fL, m = 0, 1, 2, . . . , 
when a discrete shearing-frame x-wavenumber k^m = kxf + Skyftm (or its periodic alias) 
coincides with kxf in the laboratory frame. (This discretization is the one used in a 
MHD computational code with a finite number of Fourier modes ( Yousef et al. 2008q[ ).) 
This allows the shearing-coordinate spatial average of the differenced-phase factors to 
have the requisite property for a compact ESD derivation. The resulting ESD replaces 
the time- history integrals with finite sums over m at discrete forcing times tm, and it 
replaces the continuous laboratory- frame k{t — fi) with k(t — tm)- This pathway retains 
the familiar dependence on L for a discrete Fourier series; this dependence disappears as 
L -> 00 when the shearing-periodicity pathway merges with the large-domain pathway 
as Ak and At vanish. The general behaviors of the finite-L shearing-periodicity ESD and 
L 00 ESD in (|4.22p are essentially the same. Because of the simplicity of the spatial 
averaging with the shearing-periodic boundary conditions and the analytical advantages 
of the assumpti ons of large 77 and i/, sma ll 5 7^ 0, and small tf, a. proof-of-concept ESD 
exposition is in Heinemann et all ( 2011a ). Its solution coincides with Sec. 15.1.21 Notice 
that this combined pathway achieves spatial homogeneity even without the enlarged 
ensemble of uniform mean flows in Y (Sec. 13. 5p . 
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6. General Parameter Dependences 

With the normahzation conditions p.2ip - p.22p . the non-dimensional parameters of 
the ESD equation (|4.22p are S, tj, Of, k^, 77, and 9b- A priori we are interested in 
possible dynamo behavior over their full ranges. Section lOl shows a typical "mid-range" 
example by computational integration, and Sec.[5]has analytic formulas for the parameter 
dependences of the growth rate 7 in two asymptotic limits associated with r],v ^ 00 and 

— > or 00. In this section we survey the parameters space computationally to show 
that 7 in the ESD solution is a smooth, simple function of all its parameters. 

For given parameters, a computational solution provides a particular realization of 
the random forcing in Sec. 13.31 When there is exponential growth in a fit cx e'^t 

is made over a long integration period {e.g., SAt = 10^ in Fig. [3]). The 7 value varies 
from one realization to another, but the results we report here are fairly well determined, 
as indicated by the smoothness of parameter curves based on separate estimations at 
separate parameters. Nevertheless, it is computationally laborious to obtain an ensemble 
perspective over many realizations. 

Dynamo growth occurs for finite values of < kz < kf — 1 (Fig. [S]); i.e., increasing 
kz amplifies the fluctuating helical forcing in (|4.22l) that is essential to the ESD, and 
dynamo growth is quenched by resistive decay when kz is too large. There is an optimal 
intermediate value for kz where 7 is a maximum. This behavior is approximately the 
same as evident in the analytic solutions in Sees. IS.lTD and [5.1.21 

The functional dependence of 7 on the shear S is in Fig. [71 based on optimization over 
kz with the other parameters held fixed. The dynamo growth rate increases monotonically 
with S; the slope of 7(5*) decreases for larger S. A power-law fit to 7(5*) shows an exponent 
approximately in the range 0.5-1, which is consistent with the values of 2/3 and 1 in the 
limiting formulas (|5.11|) and (|5.35l) 

The associated optimal kz{S) is always small relative to fc/ = 1, and it too increases 
with S. A power law fit shows an exponent similar to the limit values of 1/2 and 1/3 
in (|5.1ip and (|5.35l) . In the ESD there is no threshold in S for dynamo growth, given 
sufficiently small kz ^ 0. With either S — or kz = 0, there is no dynamo. 7(^2) is a 
convex function of kz that vanishes when kz is not small as well as when fc^ — >■ 0; this is 
a similar shape as in the limit formulas (|5.9p and (|5.33p . For all other parameters held 
fixed (including kz), there is a minimum threshold value of S for dynamo action, as is 
also true in the limit formulas (|5.9p and (|5.33p . 

The dependence of 7 on the forcing correlation time is in Fig. [8l again based on 
optimization over kz . 7 and kz both increase with t f . This tendency is consistent at small 
tf with the limit formulas in (j5.35p . For larger tf values the slope of 7(i/) increases with 
tf in the range surveyed here, although we know from (|5.1ip that 7 asymptotes to a finite 
value with steady forcing. The optimal kz{tf) levels off with large tf, here at a value only 
slightly smaller than kf = 1; this behavior is not anticipated by the limit formulas in 
Sec. [5| that indicate small kz for large 77. 

We demonstrate the roles of the forcing components fz and dz by alternately setting 
them to zero, dz = removes all forcing from (|4.22l) . hence has no effect on B. fz = 
retains the forcing in T but makes "H = 0; in this case B{t) shows algebraic growth 
in time but no dynamo. Thus, a dynamo requires both Uz and u± to be nonzero. By 
keeping both components non-zero but arbitrarily setting T = with H ^ in (|4.22p . 
7 is modestly increased; this confirms the interpretation of the T effect as turbulent 
resistivity that weakens dynamo growth (Sec. [5]). If J-{t) is replaced by its time-mean 
value, the dynamo behavior is essentially the same. 

Viscous and resistive diffusion both diminish dynamo growth, but they do not suppress 
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it entirely (Figs.l9lfT0|). The growth rate becomes independent of — > for fixed 77, and it 
becomes independent of 77 — > for fixed v. The latter indicates that the ESD is a "fast" 
dynamo with 7 7^ as ry — > (Roberts and Soward, 1992). At the other extreme, to 
sustain a dynamo as r/ — > 00, the value of kz{ri) must become very small so that resistive 
decay is not dominant; this is consistent with the limit formulas (jS.lip and (j5.35p . where 
7(77) decreases as a power law with exponents of -1 and -5/2, respectively. 7(1/) decreases 
with V for large v. We can take the 7/ 00 limit of (|4.22l) for general 7;, using the same 
type of approximation procedure as at the beginning of Sec. [SJ The key approximation 
in this limit is 

* dA A^l Gt (A - ii)Gl{t (6.1) 
Jo ^ 

with fi, X ^ t for the arguments of the other integrand factors. The resulting (L, z/)-limit 
mean-field equation has the same structure as (15.51) except now the electromotive force 
curl has a prefactor of I/t/ instead of 1/r/. Consequently, must decrease with large 
I' as in Fig. [HI 

The optimal 7(6*/) and kz{9f) are both largest for intermediate 9f values (Fig. [TT|) . 
The limit formulas predict a peak at df = 7r/2 and 7 = at 0/ = 0,tt {kyf — 0). 
However, these limits are based on (|5.5p after 1/ 00, which suppresses any effect of shear 
tilting in the ESD. In the more general case an up-shear orientation (0 < 0/7r/2) is more 
conducive to dynamo growth. Thus, the Orr effect of phase tilting in shearing waves (Sec. 
13. 2p augments the dynamo efficiency. This is because, when ^/ is up-shear, the helical 
forcing factor transiently increases in magnitude as (t) decreases between t — and t = 
l/tan6'//S' > when kx{t) passes through zero and thereafter becomes increasingly large 
and negative. This has the effect of transiently augmenting the effective helicity, hence 
dynamo forcing, compared to a down-shear case where |fc2;(i)| monotonically increases and 
the effective helicity only decreases with time. The magnitude of this transient dynamo 
enhancement is limited by the viscous decay that ensues during the phase tilting toward 
kx = (and beyond), consistent with the Orr effect disappearing when — 00. 

From an ensemble of numerical integrations, we find that the estimate mean value of 
7 is independent oi 9b', i-e., the initial conditions of B are not important for the dynamo 
apart from the necessity of a seed amplitude in 3 to enable the dynamo. 

The analytic solutions in Sec. l5.1l for the 77, — >■ 00 limit show that the ensemble mean 

field, £ S , has a smaller (but nonzero) dynamo growth rate 7 than the r.m.s. field for 

a steady-forcing ensemble as well as a smaller (but undetermined) 7 for rapidly-varying 
forcing. Figure [12] illustrates, for a more generic parameter set, how the components of 
the complex amplitude 13{t) vary substantially both with time and among different real- 
izations, including spontaneous sign reversals on a time scale longer than those directly 
related to the parameters (i.e., the non-dimensional fluctuation turn-over time of 1, as 
well as tf, 1/S, 1/77, and l/i^); long-interval reversals also occur for Earth's magnetic 
field. This occurs even as the mean magnetic field amplitude inexorably grows, albeit 
with evident but relatively modest low-frequency and inter-realization variability. It has 
proved to be computationally difhcult to accurately determine the ensemble mean of 13 
over many random realizations for fixed initial conditions in the general ESD (|4.22p . Our 
computational experience is consistent with the mean field magnitude typically being 
only a small fraction of the square root of the mean magnetic energy. Thus, the ESD 
with random small-scale forcing is essentially a random large-scale dynamo. 
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7. Summary and Prospects 

We derive the Elemental Shear Dynamo (ESD) model for a random barotropic force 
with a single horizontal wavevector in a steady flow with uniform shear in a large domain. 
It is a quasi-linear theory that is rigorously justified for vanishing magnetic Reynolds 
number (1/ry — > 0) and experimentally supported for more general parameters. It ro- 
bustly exhibits kinematic dynamo behavior as long as the force / has both vertical and 
horizontal components with finite forcing helicity variance; the vertical wavenumber kz 



compared to the horizontal wavenumber of the forcing; and the forcing wavenumber 
orientation is not shear-normal (i.e., kyf ^ 0). When these conditions are satisfied, the 
dynamo growth rate is larger when S is larger, the resistivity rj and viscosity v are smaller, 
the forcing correlation time i/ is larger, and the forcing wavenumber ^/ is in an upshear 
direction. The ensemble-mean of the energy of the horizontally averaged magnetic field 
grows as a dynamo, but the energy of the ensemble-mean magnetic field is much smaller. 



Reversals in ^-By (t) are common over time intervals long compared to tf. Because the 
growth-rate curves have broad maxima in both parameters and fluctuation wavenumbers 
(Sec. ini), we expect dynamo action with a broad spectrum in fc/ and kz, consistent with 
the quasi-linear superposition principle (|4.26|) . 

The ESD ingredients of small-scale velocity fluctuations and large-scale shear are 
generic across the universe, so its dynamo process is likely to be relevant to the widespread 
existence of large-scale magnetic fields. Of course, the simple spatial symmetries assumed 
in the ESD model are a strong idealization of natural flows, and the ESD is not a general 
MHD model because of its quasi-linearity assumptions. Investigation of more complex 
situations is needed to determine the realm of relevance for the ESD behavior shown 
here, especially in turbulent flows with intrinsic variability and large Reynolds number. 
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Appendix A. Derivation of Fq in (I4.17p 

This appendix fills in steps between the formal expression for the curl of the mean 
electromotive force (j4.16|) and its particular expression in the ESD (I4.17[) . Here we retain 
the convention that all vectors are horizontal. 

To provide a more compact notation, we rewrite the vertical phase factor coefficient 
for the fluctuation fleld (|4.12p as 



of the initial seed amplitude of the mean magnetic fleld ( B 




is nonzero but small 





(Al) 



where 
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+ 



- G,(t - A - m)G.(A - m) [ fc. ( - i/.(/i)S(t - A) • B{\) 

Bz X k{t — fl) 



P(A-^) 



(6,(m) (fc(A-Ai).B(A)) 



6 = 



i G„(i - Ai, A - m)G,(A - /i) [ fc, ( - irA^i)S{t - A) • B(A) 
62 X fe(/; — /i) 



P(A-/x) 



6:(m) (fc(A - /i) • S(A)) 



(A 2) 



We evaluate the three terms for ^ i^j? y in (|4.16p for each of the terms in b+ and 

b-. To do so involves the spatial average of products of factors with exponential phase 
functions, i( ± + k^z). Employing p.27p we will make use of the general identities, 



Re e'-^^")} Re {A2(m) e'-^^^^e''^^"} = ^S{p- fi) Re {^*A2e"=^^} 



Re{Ai(p)e''^('')| Re{A2(Ai)e-''^('')e"=^"| ^ ' = ^(p - ^) Re {AiAae"^^^} . 

(A3) 



The first term in ()4.16p is evaluated as follows: 



{u ■ V)6' 



Cl 
2 



= 0. 



^ dAjJI d/i [ Re{ib+(i,A,/i)e''^('^)e"=-^} 
+ Re{-i6_(i,A,M)e^'^('^^e'^^"} ] 



^0 



k^{t-p) 



(A 4) 



The formula (jASj is used to obtain the middle right-hand side, and the final result comes 
from the identify, Bz y. a ■ a — Q. 
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The second term in ()4.16p is evaluated as follows: 



{uzdz)b' 



dpG,(t-p)Re{A(p)e''^('')} fc. 



dA 



d/i 



Re 



Cl 
2 



+ Re|ib_(t,A,/^)e-''^('')e"=-^} ] J) ' 
dA J d^lG,{t - fi)k. Re { i + e"=^^ } 

d^ G,j{t - /I, A - ^i)Gu{\ - n)G^{t - fi) 

JO 

|/.P(^)fc2Re{S(t-A).B(A)e"=^-} 



Re 



k{X-fi) 



Re{iB(A)e'''-^^} 



(A 5) 



The formula (jA 3P is used to obtain the second right-hand side, and (jA 2^ is substituted 
to obtain the final result, which agrees with the first and second terms in (j4.17p . 
The final term in (|4.16p is evaluated as follows: 



{b' ■ V)u 



dA 



d^ / dpGi,{t ~ p) 



- X k{t - p) 



kit - p) 



Re 



Cl 
2 



dA 



d^G^{t - /i) 



fc2(t-p) 

{b+(t,A,M)e'*('^)e"=^-^} 

+ Re{6_(i,A,^)e~''^('^)e"^^^} ] ) Y' 
Bz X k{t — p) 



^ Rc{ii6,(p)e'*('')} 



k^{t-p) 

(^k{t~p) ■ Re{ (6*b+ +o^b_) e"^^" } 



d^Grjit - /i, A - p)Gu{X - ii)Gy{t - /i) ^''^2^Yl_ 



fc^Re|/;6^} •Rc{iS(i- A) •B(A)e'''^^}) 



(A 6) 

Again, (|A 3[) is used to obtain the second right-hand side, and (|A 2p is substituted to 
obtain the final result, which agrees with the third term in (|4.17p . In this substitution, the 
terms in b± cx 6z, 6* do not survive because they yield a factor, x k{t~ p)-k{t— p) = 0. 
This completes the derivation of (|4.17p . 



Appendix B. Computational Solution of (I4.22p 

The ESD solutions in Sec. 14.31 are obtained by numerical integration of the integro- 
differential equation system (|4.22p . This system is potentially expensive to solve because 
of the two time integrals, requiring 0{T^) operations to integrate to time T. We convert 
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this to an 0{T) system (formally comparable to the size an ODE integration, albeit 
with a much larger coefficient for T) by limiting the integration range to fixed intervals, 
t — T ^ A, ^ t and t — t ^ ii ^ X, once t > t; for smaller t values, the integrations 
start from X = fi — 0. A sufficient motivation for this approximation is that the two 
viscous decay factors (j3.19l) become vanishingly small for large values of its arguments 
t — fj, and X — II in (|4.22p . For a given v value, we determine r by the requirement that 
G'j/(t) ^ e ^ 1. In practice we typically choose e — 10~^ and make sure that the results 
do not change significantly if we further decrease the value of e. 

The domain of integration is a triangle in (A, /i)-space. Because of this, the forcing 
functions J^{fi) and only need to be retained in memory for the range t — r ^ /i ^ t 

to evaluate B{t) and advance B{t) in time. With the restricted integration intervals, the 
mean field equation (|4.22l) is 



dA 



where we have introduced the second-order matrices, 
and 



B(A), 
(Bl) 

(B2) 



Qih,t2) = -hG^{ti,t2)G,{ti)G,{t2) k-\h) + k-\t2) (e, X fc(ti)) k{t2) . (B 3) 



We may convert the double time integral in (|B ip to a double time integral 'into the past' 
via the substitutions X' = t — X and fj,' — t — ii, giving 

B{t) ^ SB^{t)ey+ dA' / dfi' T{t-fi')V{fi',fi'-X')+m{t-fi')Q{fi',f/-X') •S(t-A') • 



(B4) 

Note that, because ^ A' ^ /i' ^ r, the matrices (jB 2[) and (IB 31) can be evaluated once 
and for all in the ranges ^ t ^ t and ^ A ^ r at the beginning of the simulation. 

To discretize (IB 41) in time, we write this equation as a system of one integro-differential 
and one integral equation, viz., 



B{t) = SB4t)ey - / dXg{t,X) - Bit- X) 



and 



g{t,x)= / dfi T{t^pi)rUi,fi-x) + m{t-ii)Q{n,fi^x) 



(B5) 



(B6) 



where we have now dropped the primes from A and /i. Using the trapezoidal rule, a 
second-order accurate representation of (IB 51) is given by 



B 



n+l 



S{At) ■ B" 



At 



5(At) 



At 



m,0 



At 



4 



At 

At ^ ^„+,. 



■B 



■n+l— m 



(B7) 



where B = BijiAt) and 5"'™ = g{nAt,mAt); we have anticipated (|B86| below. The 
integer K is defined through the relation r = At. The matrix factor 5(At) arises from 
treating the shear stretching term exactly. Because (|B 7P is linear, it may be easily solved 
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for B provided the matrix 5"+^'^ can be inverted. To compute 5"'™, we again use 
the trapezoidal rule to obtain 

K-1 



_|_ ^ ^ ^■p'n—l'j;)l,l—m_^^^n—lQl,l—Tn^ 



=m+l 



0. (B86) 
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if = Re = kflm{rj. 
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Figure 1. Time series of the random forcing time variances, \ fz\^ and |/±|'^, and forcing helicity 
H{t) for a case with Of — 7r/4, 5* = 1, tf = 0.1, and v — 0.1. The discrete time step size is 
At = 0.025. 
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Figure 2. Time series for vertical and horizontal velocity variances {i.e., twice KEz(t) and 
KE±{t)) and the associated kinetic helicity response, ( u ■ lj ) (t), for random velocities 
generated by the forcings in Fig. [1] 




t 

Figure 3. Mean- field variance time series for the same case as in Figs. [Tl2l Additional case 
parameters are i] = 0.1, = 0.125, and 6b = n/4. 
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Figure 4. Steady-forcing normalized distributions D{e) for (left) |B| (solid) and the component 
magnitudes, \Bx\ (dash-dot) and \By\ (dash), and for (right) Re | Si: | (solid) and Im | By | 
(dash) for S^JeosmOf te — 6 and eo = 0.1. This case has 6 f — 8b = n/A. 




Figure 5. (Left) steady-forcing normalized distributions D(e) for \3\ at two times: 
S^/eo sinOf te = 6 (solid) and 12 (dash). Case parameters are as in Fig. |4] (Right) steady- 
forcing normalized growth rates, y^"^" and ^ as a function of evaluation time t^. 
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FlGURE 6. Mean-field growth rate 7 as a function of k^. Parameter values are S = 1, 9f — n /A, 

tf = 0.1, f = 0.1, and r; = 0.1. 
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Figure 7. Mean-field growth rate 7 associated with the optimal kz value as a function of shear 
S. Other parameters are Of = 7r/4, tf = 0.1, u — 0.1, and ri = 0.1. The dots here (and in 
subsequent figures) indicate the sampling density for this evaluation of the ESD. 
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Figure 8. Mean-field growth rate 7 associated with the optimal value as a function of forcing 
correlation time tf. Other parameters are df = 7r/4, S = 2, 1/ = 10, rj = 0.01, and At = 0.01. 
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Figure 9. Mean-field growth rate 7 associated witli the optimal fcz value as a function of 
viscosity v. Other parameters are Of = 7r/4, S = 1, tf = 0.1, r) — 0.1, and At — 0.025. 




Figure 10. Mean-field growth rate 7 associated with the optimal kz value as a function of 
resistivity rj. Other parameters are Of — n/A, S = 1, tf = 0.1, rj = 0.1, and At — 0.025. 
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Figure 11. Mean- field growth rate 7 associated with the optimal value as a function of 
forcing angle 9f. Other parameters are Of = n/4, S = 1, tf — 0.1, v — 1.0, and At = 0.025. 
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Figure 12. Five ESD realizations in different colors of the mean field component vari- 
ances, \Ba^f{t) and \Byf{t) (top row), and of their real and imaginary parts normalized by 
Brma = (bottom two rows). Parameters are kz — 0.14, S — 1, Of = it/A, tf = 0.1, u = 0.1, 

and 77 = 0.1. All realizations have the same initial condition B(0). 



